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Abstract. We discuss the internal structure of stars in the mass range 1.5 to 4Mq from the 
PMS to the subgiant phase with a particular emphasis on the convective core and the convective 
superficial layers. Different physical aspects are considered such as overshooting, treatment of 
convection, microscopic diffusion and rotation. Their influence on the internal structure and on 
the photospheric chemical abundances is briefly described. 

The role of binarity in determining the observed properties and as a tool to constrain the 
internal structure is also introduced and the current limits of theories of orbital evolution and 
of available binary data-sets are discussed. 
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1. Introduction 

The theoretical evolution of A stars is extremely simple to compute if one ignores 
complex phenomena such as gravitational settling, radiative forces, rotation, turbulent 
mixing, magnetic fields, binarity... We shall briefly introduce the infuence of different 
treatments of convection, diffusion and rotation after a presentation of a "conservative" 
situation in which the only mixings come from convection and overshooting. Binarity as 
a tool to constrain stellar models is then discussed. 

2. A type stars: evolution 

As A stars are located on or near the main sequence, we shall start this discussion 
with the core hydrogen burning phase. 

2.1. Core hydrogen burning 

Figuren(lp) shows the evolutionary track of a 3M© with a solar chemical composition. 
The evolution with time of the radiative gradient, given by 

Vrad ~ — K (2.1) 

m 

is shown in figure [21 (lp) in the inner 30 % of the mass for selected models on the 
evolutionary path. 

The adiabatic gradient, Vad, is also drawn on this figure and its decrease towards the 
centre shows that the influence of the radiation pressure increases with time. Starting 
from the right of the figure, the radiative gradient increases when reaching the CNO 
burning central layers, due to the increasing value of L/m, and hydrogen burning thus 
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Figure 1. Left panel (lp): Evolutionary track for a 3Mq star computed without overshooting. 
Right panel (rp): Same figure with an overshooting parameter of 0.2 showing a wider MS track 
than in the left panel 




m/M m/M 

Figure 2. Left panel (lp): Radiative gradient for different models on the main sequence (black 
curves - first two dots starting from the left in figure and in the early H shell burning phase 
(gray curves - last two dots starting from the left in figure 0. Right panel (rp): Same figure 
with an overshooting parameter of 0.2 

takes place in a convective core. Inside this core, the hydrogen abundance, X, decreases. 
The opacity, k, being proportional to (1+X), decreases as time goes on. This means that 
the convective core has its maximum extension in mass on the ZAMS and then slowly 
recedes as hydrogen is transformed into helium. The difference in X between the core 
and the non-burning regions thus increases with time without any discontinuity at the 
edge since the convective core recedes. A gradient of hydrogen builds up and a plateau 
in Vrad appears (see figure |21 (lp)). 

The decrease of X leads to an increase in the temperature (Figure |3| (lp)) which is 
limited to the central layers while the other layers become cooler and cooler. The central 
part of the star is contracting while the other layers are expanding. Once X comes close 
to zero in the core, the increase in temperature becomes more pronounced (see figure [21 
(lp)). This is done by a global contraction (the so-called second gravitational contraction). 

2.2. Hydrogen shell burning 

In the fi- gradient zone, a hydrogen burning shell starts to develop. A secondary maximum 
is seen in Vrad (Figure|3](lp)) as well as in the nuclear enregy rate, e (Figure0](lp)). Soon 
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Figure 3. Left panel (lp): Temperature versus mass for different models on the main sequence 
(black curves - first two dots starting from the left in figure Q and in the early H shell burning 
phase (gray curves - last two dots starting from the left in figure 0. Right panel (rp): Same 
figure with an overshooting parameter of 0.2 
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Figure 4. Left panel (lp): Nuclear energy rate versus mass for different models on the main 
sequence (black curves - first two dots starting from the left in figure and in the early H shell 
burning phase (gray curves - last two dots starting from the left in figure 0. Right panel (rp): 
Same figure with an overshooting parameter of 0.2 



after, Vrad becomes smaller than y a d and the convective core vanishes. As the energy 
production stops, an isothermal core appears. The models are now located between the 
last two dots on the evolutionary track in figure ^ (lp) . As hydrogen is burned in the 
shell, the isothermal core increases in mass. However, its mass cannot exceed the limiting 
mass of Schonberg-Chandrasekhar, given by 

(£) ^0.37^ (2.2) 
VAf/sc m 

which is of the order of 0.08. When this value is reached, the core must contract in 
order to create a temperature gradient (see figure (lp)) and this is accompanied by an 
expansion of the envelope where the temperature decreases. 

2.3. Effect of overshooting 

In order to fit CM diagrams of open clusters as well as eclipsing binaries, overshooting 
seems to be needed with an extent increasing with stellar mass. We have computed a 
similar evolution with an overshooting parameter of 0.2. A larger hydrogen reservoir 
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Fig ure 5. Left pane l (lp): P MS tracks (see text) wi th two birthlines (bottom one: 
|Palla fc Stahler (1993)] top one: |Behrend fc Maeder (2001)1 Center panel (cp): Evolutionary 
tracks computed with FST (dashed curves) and with MLT (a = 1.6) (faded curves). Right 
panel (rp): Fraction of energy transported by radiation versus logT in the external layers. (1) 
With FST : Only one (the smallest) minimum near logT = 4.1 (Hi and He i) and no mini- 
mum near logT = 4.7. (2) With MLT (a = 0.5) : One minimum (the intermediate one) near 
logT = 4.1 and no minimum near logT = 4.7. (3) With MLT (a = 1.75) : One (the largest) 
minimum near logT — 4.1 and a second small one near logT = 4.7 (He n) 



due to the increase of the mixed region evidently translates into a longer core hydrogen 
burning phase as well as in a longer track in the HR diagram (figure ^ ( r P))- However, 
after the turn-off, for a similar variation in effective temperature, the evolution with 
overshooting is much more rapid. Figure [21 (rp) shows that Vrad never reaches zero 
so there is no plateau in temperature (figure |21 (rp)). The reason lies in the fact that 
with the increase of the mixed region, the exhaustion of hydrogen takes place in a core 
whose mass is already higher than (jj) sc - This results in a much quicker contraction 
of the central regions. The star moves more rapidly and for a similar change in effective 
temperature, the shell, shown by a maximum in the distribution of e, is still at the same 
mass fraction with overshooting while the shell has already moved significantly in mass 
without overshooting (compare lp with rp in figure 0J. 

Although the ^.-gradient has a smaller slope with overshooting than without overshoot- 
ing, the presence or absence of an isothermal core leads to a similar extent of the helium 
core at a similar effective temperature. 



2.4. Pre-main sequence evolution 

Let us now turn toward the pre-main sequence evolution (figure|5j). The two quasi- vertical 
lines on the Hayashi tracks and soon after show where the star becomes partly radiative 
and then completely radiative. The loops near the MS are the signatures of the CNO reac- 
tions evolving toward equilibrium. We have checked the effect of the new (smaller) cross- 



section ofthe slowest ofthe CNO reactions, i.e. 14 N(p,7) 15 jFormicola fc LUNA (2004) \ 



Although the approach toward equilibrium is slower, the difference in the main sequence 
lifetime is not significant. 

The points in the HR diagram where the stars become visible define the birthline. Two 
birthlines are drawn in figure|Sl(lp)). The bottom one comes from Palla & Stahler (1993) 
and is constructed assuming a constant accretion rate of 10~ 5 M Q /yr. This value was 
obtained by fitting the Herbig Ae and Be objects in the HR diagram. The top one comes 



from an analysis by Behrend & Maeder (2001) where the accretion rate 
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disk 



(2.3) 
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Figure 6. Illustration of diffusion. See text 



has been derived in order to reproduce the observations in the HR diagram, in particular 
for the most massive stars. The time spent on the Hayashi track depends on the adopted 
birthline. This results in a different transfer of angular momentum from the star to the 
disk, which in turn can change the angular momentum at the beginning of the MS phase. 

In these fully or partly convective phases, the treatment of convection affects the 
location in the HR diagram, creating a significant difference in effective temperature. 
Figure (cp) shows the location of the Hayashi tracks for models computed with the 
FST treatment of convection (dashed curves) and with the MLT treatment (a = 1.6) 
(faded curves). The difference in effective temperature is of the order of 200 K. As the 
convective envelope recedes, the tracks become undistinguishable. 

However, the extremely thin surface convection zone which remains on the main se- 
quence is affected (figure 03 (rp)). The most superadiabatic the temperature gradient 
is, the less efficient is the convection. For a 1.8 Mq, it can make the convective He n 
ionization zone appear or disappear. This is important since the thickness of the mixed 
superficial layers is crucial in explaining the abundance anomalies. 

2.5. Gravitational settling 

With such thin convective envelopes, gravitational settling is very efficient. Figure El 
illustrates this point. 

In the left panel, a duck-star has all its chemical elements (chicken) showing at the 
surface. In the middle panel, the MS phase is symbolized by the crossing of a grid. 
If radiative forces are not taken into account together with gravitational settling, the 
duck-star ends its MS phase with only one element (chick) remaining at the surface. 

However, radiative forces can in turn be too efficient and another mechanism must enter 
the game. It is the turbulent mixing which will not only affect the surface abundances 
but will induce changes in the internal structure as well, especially in the thickness 
of the convective envelope ( |Richard, Michaud &: Richer, 2001) 1 . One striking example is 
the formation of a convective iron zone at a temperature of 200000 K caused by an 
accumulation of iron-peak elements as a result of radiative forces. According to the mass 
of the star, this zone can merge with the H-He ionisation zone during the main sequence 
or remain detached, which affects the surface abundances of the iron-peak elements. 

3. A type stars: rotation 

From their location in the HR diagram, A-type stars are expected to rotate rapidly and 
not to be affected by magnetic braking. They can reach indeed rotational velocities up to 
300 kms -1 . Abt & Morrel (1995, hereafter AM95) showed however, that the distribution 
of rotational velocities has a bimodal shape, with virtually all the Am and Ap stars having 
equatorial rotational velocities (w ro t) less than 120 kms -1 , and most of the normal AO- 
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FO main-sequence stars having v ro t > 120 kms 1 . They concluded (see also Abt (2000)] ) 



that rotation alone can explain the occurrence of abnormal or normal main-sequence A 
stars, and that the apparent overlap between their u ro t distributions is only due to our 
inability to distinguish marginal Am stars from normal ones, or to disentangle rotational 
and evolutionary effects. 

A large number of observational data indicates that during MS evolution of A-type 
stars: 1) mass loss is limited to ~ 2.10~ 10 M Q yr _1 . (Lanz & Catala, 1992); 2) rota- 
tional velocity does not depend on age; and 3) no significant angular momentum loss 
by magnetic braking is observed (e.g. Wolff & Simon , 1997; Hubrig, North & Medici, 
2000). Consequently, the observed v ro t distribution must be determined by the angular 
momentum evolution during the pre-main sequence (PMS) phase. In a study of the dis- 



tribution of angular momentum in Orion stars, Wolff, Strom & Hillenbrand (2004) find 
this hypothesis consistent with a simplified model of PMS evolution in which angular mo- 
mentum is lost by interaction with the protostellar disk during the convective phase, and 
conserved in the radiative one. Furthermore a core-envelope decoupling occurs during the 
convective-radiative transition (see also |St§pieri, 200 and StQpie n~&; Landstreet, 2002| for 
Ap PMS angular momentum evolution). The understanding of angular momentum evo- 
lution during the PMS and the MS is fundamental to understand the Ap and the Am 
phenomena. 

3.1. Modeling the evolution of a rotating star 

Rotation has several different effects on stellar evolution: change of the internal hydro- 
static equilibrium, changes in the apparent effective temperature and luminosity (Perez 
Hernandez et al. 1999 and references therein); transport of chemicals and of angular 
momentum by shears in differentially rotating stars, by meridional circulation an by 
horizontal turbulence (e.g., Zahn, 1974; Knobloch & Spruit, 1982). 

The effect of rotation on the stellar evolution has been treated with different ap- 
proaches: Endal & Sofia (1981) assigned to each transport process (including meridional 



circulation) a diffusion coefficient, while [Zahn (1992a)] and Maeder & Zahn (1998) made 
the hypothesis that differential rotation in the radiative zone of a non-magnetic star 
gives rise to anisotropic turbulence (much stronger in the horizontal direction than in 
the vertical one due to the stratification) , and those result in "shellular" rotation. In this 
model the effective diffusion coefficient for the chemicals is: 



\rU(r)\ 2 , r 1 d(pr 2 U) U d\nr 2 fl 



(3.1) 



3pr dr 2 dlnr 

Where U is the meridional circulation velocity, Ch is a free parameter related to the turbu- 
lent horizontal viscosity, and the other quantities have the habitual meaning. |Palacios, Talon, Charbonnel, et al. ( 2003) 
has applied the Zahn's modeling of rotation with the formulation by |Maeder fc Zahn (1998)] ~ 
(including also microscopic diffusion but not radiative accelerations) to masses from 1.35 
to 2.2 Mq. We recall here two results of this application: 1) the rotation profile inside 
the star shows differential rotation mainly close to the convective core. 2) The models 
(computed without overshooting), show a wider MS when rotation is included. However, 
while the lower mass models increase their MS lifetime by a 20%, the highest masses 
increase it by only 10% (see their Fig. 1). This trend is opposite to that observed in open 
clusters and in binary systems, where the fitting of observations requires to increase 
the overshooting parameter with the mass (i.e. Andersen, Nordstrom, &: Clausen (1990)| 
Ribas, Jordi & Gimenez, 2000). 



Recently, Maeder (2003) and independently Richard & Zahn (1999) have updated the 
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value for the horizontal diffusion coefficient, Ch by a factor of 10 2 with respect to the value 
used in Maeder & Zahn (1998) Its implementation in a stellar evolution code (Mathis, 
Palacios & Zahn 2004) leads to: 1) enhanced mixing, and 2) significant changes in the 
profile of chemical mixing with depth (see their Fig. 2). To which degree will the new D e g 
(with a non-zero value until the boundary of the convective core) affect the MS width 
and the surface chemical abundances? 

The scheme proposed by AM95 matches qualitatively well the predictions of the mi- 
croscopic diffusion models (e.g. Richer, Richard & Michaud 2000), in the hypothesis that 
the extra-mixing required to decrease the microscopic diffusion is induced by rotation. 
Some problems are however left, for instance: 1) the "normal" late B and early A-type 
main-sequence stars show a true bimodal v mt distribution (Royer, Gomez & Zorec, these 
proceedings); 2) the binary V392 Carinae has a usini of 27 kms" 1 and no peculiar 
abundance (Debernardi & North 2001). 3) There is no correlation between the strength 
of chemical peculiarities and i> rot (Ersp amer fc North, 2003| |. Models including micro- 
scopic diffusion and radiative accelerations show that the chemical abundances are very 
sensitive to the thickness of the mixed layer (e.g., Alecian 1996, Hui-Bon-Hoa, 2000 
|Richard, Michaud fc Richer, 2001| ), so that we should expect a signature of t> rot on the 
abundances, if rotation is the responsible of the extra-mixing below the convective enve- 
lope. A first approach to explicitly include a model of rotationally induced mixing and a 
complete treatment of microscopic diffusion is in progress (Richard, Talon & Michaud, 
these proceedings). 

Finally, there are also other effects that should be taken into account in the study of 
the evolution of a rotating A-type star: 

1. Interaction convection-rotation. Rotation can reduce the efficiency of convection and 
is able to reduce the extension of the overshooting region (e.g. |Julien, Werne, Legg et al. 1997| ). 
Recent numerical simulations of the a rotating convective core (Browning, Brun & 
Toomrc, 2004) show that rotation leads to a variation of convective overshooting and 
penetration, and induces internal waves, meridional circulation, and differential rotation 
at the core boundary. 



2. Magnetic field. Maeder & Meynet (2003) have shown that the Tayler-Spruit mag- 
netic instability (Spruit 2002) could take place in the interior of stars with small magnetic 
fields and differential rotation, induce a process of angular momentum transport much 
more efficient than that due to the meridional circulation and horizontal turbulence and 
lead to solid body rotation on a short timescale. Furthermore, in Ap stars, which have 
suffered strong braking during their PMS phase (St§pieh, 2000), differential rotation 
together with magnetic field could lead to magneto-rotational instabilities. These could 
transport angular momentum from the interior to the surface ( Arlt, Hollerbach, fc Rudiger (2003)) 



and Artl, these proceedings), and result in solid body rotation of these stars, as expected 
from some observational data (Hubrig, North & Medici, 2000). 



4. A type stars: binarity 

Binarity plays a fundamental role in the origin and definition of A-type stars chem- 
ical peculiarities. The pioneering work of Abt (1961, 1965) showed that most (he even 
suggested all) Am type stars are relatively short period binaries and that the period 
distribution of binaries with Am and with non-peculiar A-type components are com- 
plementary. Normal A-type components are found in systems with period shorter than 
~ 2.5 or longer than ~ 100 days, while Am stars form binaries with period in the range 
2.5-100 days. 

Abt's statements have somewhat lost of strength with time, and with the increas- 
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ing size of studied binary samples. Some overlap between the period distribution of 
normal and Am binaries was found by the same author (AM95) and the binary fre- 
quency among Am stars has steadily decreased with time (see Abt & Levy (1985) and 
North, Ginestet, Carquillat, et al. (1998)||Debernardi (2000)| for results based on CORAVEL 



surveys) . The latter author finds a frequency as low as 57% (but CORAVEL samples are 
certainly biased against fast rotation and, because of the limited time span of the survey, 
against longer period-eccentric orbit binaries). The - at any rate - high frequency has 
a straightforward explanation: in close binaries the spin-orbit synchronization by tidal 
mechanisms can efficiently brake the stellar rotation to values compatible with the Am 
phenomenon. 

The other relevant connection between binarity and A-type star peculiarity is in the low 
frequency of spectroscopic binaries among magnetic Ap stars ( Abt fc Snowden (1973)[ 



Gerbaldi, Floquet, fc Hauck (1985)[|North, Ginestet, Carquillat, et al. (1998)| ). The cur 



rent explanation is that the strong magnetic fields prevent close binary formation (how- 
ever 



Budaj (1999) suggested instead that it is binarity to affect magnetisms and not the 
other way around). 

While the general outline of the binarity-peculiarity connection is well established, 
there are still shortcomings in the theory of tidal synchronization. Besides, possible draw- 
backs in the interpretation of the observations are caused by the origin and composition 
of the observed samples (biases, selection effects). 

Two competing theories were developed to explain the observed levels of orbital circu- 



larization and spin-orbit synchronization in close binaries: the tidal theory of Zahn (1992b, and references therein) 



and the hydrodynamical theory of Tassoul & Tassoul (1992, and references therein) The 



necessary ingredients of the first are tidal bulges and an efficient dissipation mechanisms; 
in absence of synchronism, dissipation causes a lag of the bulge and hence a torque, which 
tends to establish synchronization (and orbit circularization) . The dissipation mecha- 
nisms at work are different for late and early type stars, in the first case it is turbulent 
dissipation in the convective envelope retarding the equilibrium tide, in the second it is 
radiative damping acting on the dynamical tide (forced gravity waves are emitted from 
the lagging convective core and are damped in the outer layers). 

Tassoul's theory is based, instead, on the idea that while the torque due to dissipation 
processes is negligible, transient strong meridional currents are produced by the tidal 
action and transfer angular momentum between the stellar interior and an "Ekman 
layer" close to the surface. As a consequence if the rotation period is shorter than the 
orbital one, the star is braked. 

Both theories yield time-scales for synchronization (t ms ) and circularization (t mc ). For 
early type stars, in Zahn's case: 

1 / GM \ ^ 2/ -,MR 2 ( 1 ( GM \ ^ , ,n f R\^ 

(4.1) 

where q and a are the mass ratio and the semi-axis of the binary, M, R, I the mass, the 
radius and the inertial moment of the braked star and E2 a constant related to the size 
of its convective core. In Tassoul's theory: 



1 iv . i/ML 2 \Yfi\ 8 1 « ^, , 11 _o / ML 2 \ */ 

(4.2) 

with L being the star luminosity, /? its fractional gyration radius, N the ratio between 
eddy and radiative viscosity (N — for radiative envelopes). The other factor, 7, was 
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Figure 7. Left panel: the rotational velocities (v sin i) of spectroscopic binaries with Am com- 
ponents versus orbital period. The hyperbolic curves give the theoretical relation (i = 90°) 
for synchronized systems of 2 Mq and different radii. The di scontinuous dotted line is the lower 
boundary of the "avoidance zone" , according to |Budaj (1996) 1 the vertica l lines are the expected 
upper boundary for circularization, according to llNortnTT^ahn, 2003] dotted, R = 2.1 Rq, 
3=1; dash-dotted, R = 3i?Q, q = 0.2. Right panel: the longest orbital period of a SB1 binary, 
corresponding to the assumed minimum observable radial velocity amplitude K\ = lOkrns - , 
is plotted as function of mass ratio and eccentricity for primary mass M = 2 Mq and inclination 
i = 90° 



introduced by Claret, Gimenez, & Cunha (1995) to take somehow into account the fact 
that the orbital evolution is derived (as in Zahn's case) under the hypothesis of small 
deviations from synchronism and from circular orbit. The circularization time is in general 
two/three orders of magnitude longer than that of synchronization, due to the larger 
amount of angular momentum stored in the orbit. It is evident by comparing Eqs. 14.11 
with 14.21 that Tassoul's mechanism has a longer range and a much higher efficiency, 
especially for early-type stars. The timescales are - at any rate - only an indication 
of the process speed, and cannot replace the full integration of the orbital evolution 
equations. 

In principle the comparison between the expected and the observed degree of syn- 
chronization/circularization in binaries of known age and accurate dimensions could 
discriminate between the two theories, and early type stars are the best choice for 
such a test. However, t he detailed treatment of Claret, Gimenez, fc Cunha (1995) and 
Claret & Cunha (1997) (with simultaneous computation of orbital and stellar evolution 
and a comparison sample of observed binaries with the best known parameters) showed 
that the results are still inconclusive, as both models can explain only a part of the ob- 
servations. Zahn's mechanism is not efficient enough for early type stars (but the pre-MS 
phases were not considered in the abovementioned papers) while Tassoul's one is too 
efficient, unless a high value of 7 is introduced (to shorten the time-scales by a factor 
~ 40). The most promising scenario remains therefore that suggested in a study of circu- 
larization in binaries with A-type components by [Matthews &: Mathieu (1992) a process 



composed by an important PMS phase of orbital evolution followed by a MS phase, which 
however still awaits full modeling. 



The recent survey of A-type star rotational velocities by Royer et al.(2002), hereafter R02 
provides excellent material for the study of synchronization in binaries with Am com- 
ponents. Figure {7\ shows the updated vsmi from R02 versus orbital period, for the Am 
binaries listed in |Budaj (1996)[ and the relations expected in case of spin-orbit syn- 
chronization. The vertical lines are the periods corresponding to the fractional radius 
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(R/a) that should assure circularization ( |North fc Zahn, 2 003) in the case of a 2M Q 
MS-primary. The data are generally in good agreement with the expectations, with a 
marginal indication of radii larger than in normal stars (as the synchronization relations 
should be upper envelopes, having assumed i = 90°). 

The R02 velocities are generally higher than those measured by AM95 and used by 
Budaj; the reason is a systematic effect in the velocity of standard stars used by AM95. 
As a consequence the "avoidance zone", found by |B uda j~ 996) for periods 4 < P < 20 



days (i.e the region above the dotted line in Fig. 0) is no longer empty and it is not 
necessary an ad-hoc mechanisms to explain it, such as the "tidal mixing" introduced by 
the abovementioned author; the decreasing number of systems found in proximity of the 
synchronism lines can be due to a dependence on period of the braking efficiency. 

The other feature appearing in Figure the lack of systems between 180 < P < 800 
days, is quite probably a selection effect due to the low probability of discovery as a 
spectroscopic binary. An idea of the relevance of selection effects can be derived from 
the right panel of Fig. This shows, for a M = 2 M primary, the largest orbital 
period corresponding to a given minimum observable radial velocity amplitude, K\. This 
is a function of mass ratio, q, and eccentricity, e (in the hypothesis of inclination i — 
90°). Even in the conservative hypothesis K\ = 10 kms^ 1 (but Hog eveen (1992)| on the 



basis of a detailed study of selection effects for spectroscopic binaries, suggests 25 kms -1 
for A-type stars) Figure indicates that the largest period of binaries detectable by 
spectroscopic surveys is of some hundred days (and many systems will escape detection, 
as the distribution of mass ratio of single lined spectroscopic binaries is peaked around 
q = 0.2). The gap in the period distribution could, therefore, be due to the transition 
between spectroscopic binaries, discovered as such, and visual binaries with known radial 
velocity. 

In conclusion, a proper treatment of selection effects is a main issue to solve if we want 
to extract correctly the rich information that A-type stars can provide on binary secular 
evolution and that, in their turn, binaries can yield on A-type star properties. 
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